Stability evaluation method and system of direct-drive wind turbine generator

ABSTRACT

The application involves a stability evaluation method and system of direct-drive wind turbine generator, which belongs to the technical field of wind power generation. It solves the problems of poor stability of direct-drive wind turbine generator, low safety performance, failure to realize mutual cooperation between online stability evaluation and parameter adjustment in the existing technology. Firstly, the variation of direct-drive wind turbine generator terminal voltage, current, power and PLL angle of direct-drive wind turbine generator is measured, the terminal energy and the energy negative gradient of direct-drive wind turbine generator is calculated. Then, the system stability is assessed according to the value of energy negative gradient, the influence of the critical parameters on system stability is analyzed and the preliminary adjustment strategy is proposed. It realizes the mutual cooperation between stability level evaluation and parameter adjustment, and improves the working efficiency of the system.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority to Chinese Patent Application No. 201910695024.6, filed on Jul. 30, 2019, which is hereby incorporated by reference in its entirety.

TECHNICAL FIELD

The application relates to the technical field of wind power generation system, and particularly relates to a stability evaluation method and system of direct-drive wind turbine generator.

BACKGROUND

As a renewable, clean, pollution-free and inexhaustible energy, wind energy has become an important direction of the development of renewable energy in the world. The direct-drive wind turbine generator has become the key models of research because of its advantages of small mechanical loss, high operation efficiency and low maintenance cost and the like.

The problem of sub/super-synchronous oscillation caused by the direct-drive wind turbine generators connected to weak power grid has gradually become prominent, and has occurred even more frequent especially in collecting areas of wind farms. Therefore, it has become an urgent problem to study how to assess the stable level of system fast and effectively in the early stage of oscillation and provide accurate theoretical reference for oscillation suppression and system dispatch.

At present, the stability analysis of sub/super-synchronous oscillations has gradually become the focus of many experts and scholars at home and abroad. However, existing analysis methods have a large calculation workload for solving nonlinear differential equations, they can only provide offline stability criteria, it is difficult to perform real-time online stability assessment, and it is impossible to perform effective online parameter adjustment.

SUMMARY

Based upon former analysis, the present embodiment aims to provide a stability evaluation method and system of direct-drive wind turbine generator to solve the problems of poor stability, low safety performance, and failure to realize online stability evaluation and parameter adjustment in the existing technology.

The aims of the application are mainly realized by the following technical scheme.

A stability evaluation method of direct-drive wind turbine generator includes the following steps:

measuring the variation of voltage, current, active power and PLL angle at the terminal of direct-drive wind turbine generator in one oscillation period

$\frac{2\;\pi}{\omega};$

based on the measurements of the variation, calculating the terminal energy of direct-drive wind turbine generator with the established terminal energy model of direct-drive wind turbine generator;

according to the calculated terminal energy of direct-drive wind turbine generator, obtaining energy negative gradient of direct-drive wind turbine generator, and assessing system stability level according to the value of the energy negative gradient; and

based on the influence of critical parameters of stable operation points of direct-drive wind turbine generator, PLL and transmission lines on the energy negative gradient, obtaining the adjustment law of critical parameters affecting the stability of the system.

On the basis of the above scheme, the application also makes the following improvements:

furthermore, the terminal energy model of direct-drive wind turbine generator is:

${\Delta{W_{PMSG}(t)}} = {{\int_{t}^{t + \frac{2\pi}{\omega}}{\Delta P_{e}d{\Delta\theta}_{pll}}} + {\int_{t}^{t + \frac{2\pi}{\omega}}\left( {{\Delta i_{d}d\Delta u_{q}} - {\Delta i_{q}d\Delta u_{d}}} \right)}}$

Where ΔP_(e) is the variation of direct-drive wind turbine generator output active power compared with the steady-state value, Δθ_(pll) is the variation of PLL angle, Δi_(d), Δi_(q), Δu_(q) and Δu_(d) are respectively the variation of d-axis and q-axis components of current and voltage at the terminal of direct-drive wind turbine generator, ω is the oscillation frequency of dq axis.

Furthermore, the calculation process of energy negative gradient includes:

step 1: discretizing the terminal energy:

${\Delta{W_{PMSG}(i)}} = {\sum\limits_{t = t_{i}}^{t_{i} + \frac{2\pi}{\omega}}{\Delta{{W_{PMSG}(t)}.}}}$

The value of terminal energy of direct-drive wind turbine generator at the time t₁, . . . t_(n), i.e. ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) can be obtained; n is the number of data.

Step 2: calculating the energy negative gradient η with discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)).

Furthermore, the step 2 includes:

calculating the discretized average value of the discretized terminal energy of direct-drive wind turbine generator by the following formula;

${\Delta{\overset{\_}{W}}_{PMSG}} = {\frac{1}{n}{\sum\limits_{i = 1}^{n}{\Delta{W_{PMSG}(i)}}}}$

with the discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) and its average value, obtaining the energy negative gradient by the following formula:

$\eta = {- {{\nabla\left\lbrack \frac{❘{{\Delta{W_{PMSG}(t)}} - {\Delta{\overset{\_}{W}}_{PMSG}}}❘}{\left( {\frac{1}{n}{\sum\limits_{t = 1}^{n}\left( {{\Delta{W_{PMSG}(t)}} - {\Delta{\overset{\_}{W}}_{PMSG}}} \right)^{2}}} \right)^{\frac{1}{2}}} \right\rbrack}.}}$

Where ∇[⋅] is the symbol of gradient.

Furthermore, the stability level of the system can be assessed according to η: When η is positive, the system is stable and the bigger the value of η is, the higher the system stability level is; When η is zero, the system is critically stable; When η is negative, the system will go unstable, and as the absolute value of η increases, system stability level will drop.

Furthermore, the adjustment law of critical parameters affecting the stability of system can be obtained through the following steps:

step S1: applying the PLL dynamic function to the energy negative gradient expression, and obtaining the energy negative gradient with direct-drive wind turbine generator stability operation points and PLL critical parameters:

$\eta_{2} = {{- \left( {{I_{{sd}0}\omega\frac{\cos\varphi_{c}}{2}} + {I_{sq0}\omega\frac{\sin\varphi_{c}}{2}}} \right)}{\text{⁠}{\left\lbrack {{{A_{-}(t)}\Omega_{u -}\alpha_{u}e^{\alpha_{u}t}} + {\frac{d{A_{-}(t)}}{dt}\Omega_{u -}e^{\alpha_{u}t}}} \right\rbrack|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}\text{⁠}{- {U_{{sd}0}\omega{\frac{\cos\left( {\varphi_{c} - \varphi_{RL}} \right)}{2}\left\lbrack \text{⁠}{{{A_{-}(t)}\alpha_{i}\Omega_{i -}e^{\alpha_{i}t}} + {\frac{d{A_{-}(t)}}{dt}\Omega_{i -}e^{\alpha_{i}t}}} \right\rbrack}}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}\text{⁠}{{+ \left( {{I_{{sd}0}\omega\frac{\cos\varphi_{c}}{2}} + {I_{sq0}\omega\frac{\sin\varphi_{c}}{2}}} \right)}{{\left\lbrack {{{A_{+}(t)}\Omega_{u +}\alpha_{u}e^{\alpha_{u}t}} + {\frac{d{A_{+}(t)}}{dt}\Omega_{u +}e^{\alpha_{u}t}}} \right\rbrack|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}\text{⁠}{{+ U_{{sd}0}}\omega{\frac{\cos\left( {\varphi_{c} - \varphi_{RL}} \right)}{2}\left\lbrack \text{⁠}{{{A_{+}(t)}\alpha_{i}\Omega_{i +}e^{\alpha_{i}t}} + {\frac{d{A_{+}(t)}}{dt}\Omega_{i +}e^{\alpha_{i}t}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}\text{⁠}{{- 2}\omega K_{pll}^{2}\Omega_{u -}\Omega_{u +}e^{\alpha_{u}t}\alpha_{u}\sin{\varphi_{c}\left\lbrack {{I_{{sd}0}{\sin\left( {\varphi_{u +} + \varphi_{u -}} \right)}} - {I_{sq0}{\cos\left( {\varphi_{u +} + \varphi_{u -}} \right)}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}.}}}}}}$

Where A⁻, A₊ respectively represent the amplitude of PLL angle variation Δθ_(pll) in sub-synchronous frequency and super-synchronous frequency; I_(sd0), I_(sq0), U_(sd0), U_(sq0) are the original value of d-axis and q-axis components of current and voltage at the terminal of direct-drive wind turbine generator; α_(u) and α_(i) are the attenuation coefficients of oscillation components of voltage and current; Ω_(u−), Ω_(u+), Ω_(i−) and Ω_(i+) are respectively the amplitude coefficients of sub/super synchronous components of voltage and current; φ_(u−) and φ_(u+) are the initial phase angles of sub-synchronous and super-synchronous voltage; φ_(i−) and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents; Ω_(u+)=Ω_(i+)X₊, Ω_(u−)=Ω_(i−)X⁻, φ_(u−)=φ_(i−)+φ_(RL) and φ_(u+)=φ_(i+)+φ_(RL), where φ_(RL) is the phase angle of transmission lines; k_(pθ) and k_(iθ) are the proportion gain and integration gain of PLL;

${\varphi_{c} = {\arctan\frac{\omega k_{p\theta}}{k_{i\theta}}}},{{K_{pll} = {❘\frac{{k_{p\theta}s} + k_{i\theta}}{s^{2} - {U_{sq0}\left( {{k_{p\theta}s} + k_{i\theta}} \right)}}❘}_{s = {j\omega}}};}$

Applying the dynamic function of system transmission lines to the energy negative gradient expression, and obtaining the following energy negative gradient with critical parameters of transmission lines:

$\eta_{1} = {{{- 2}\Omega_{i +}\Omega_{i -}e^{{2\alpha},t}\alpha_{i}X_{L}{{\omega sin}\left( {\varphi_{i +} + \varphi_{i -}} \right)}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}\text{⁠}{{+ {}}{{{2\left\lbrack \text{⁠}{{\Omega_{i -}^{2}e^{{2\alpha},t}\alpha_{i}\omega R} - {\Omega_{i +}\Omega_{i -}e^{{2\alpha},t}\alpha_{i}R{{\omega cos}\left( {\varphi_{i +} + \varphi_{i -}} \right)}}} \right\rbrack}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}.}}}}$

Where Ω_(i−) and Ω_(i+) are the amplitude coefficients of current sub-synchronous component and super-synchronous component respectively; α_(i) is the attenuation coefficients of oscillation components of current; R and X_(L) are line resistance and the equivalent reactance of the transmission lines; φ_(i−) and φ_(i+) are respectively the initial phase angles of sub-synchronous and super-synchronous currents;

Step S2: based on the constructed energy negative gradient with the critical parameters of direct-drive wind turbine generator stability operation points and PLL, and the one with the transmission lines critical parameters, obtaining the influence pattern of these parameters on the energy negative gradient;

The influence patterns include: reducing I_(sd0) or increasing I_(sq0) helps to improve the stability level of system; The bigger R is and the smaller XL is, the more stable the system is; Reducing PLL parameters Kp and Ki, the system stability can be enhanced.

On another hand, the present embodiment proposes a stability evaluation system of direct-drive wind turbine generator, which includes a PMU, a first processor, a second processor, a third processor and result output terminal;

The PMU is used to collect the variation of voltage, current, active power and PLL angle at the terminal of direct-drive wind turbine generator;

The first processor is used to calculate the terminal energy of direct-drive wind turbine generator;

The second processor is used to assess the stability level of direct-drive wind turbine generator according to the value of energy negative gradient;

The third processor is used to analyze the influence pattern of direct-drive wind turbine generator stable operation points, PLL and transmission lines critical parameters on the energy negative gradient of direct-drive wind turbine generator and propose the adjustment law of critical parameters to improve the stability of the system;

The result output terminal is used to output the stability level of the system and adjustment law of the critical parameters affecting system stability.

Furthermore, the PMU equipped at the terminal of direct-drive wind turbine generator and PLL collects the variation of voltage, current, active power and PLL angle.

Furthermore, the first processor receives the data transferred by the PMU, the terminal energy of direct-drive wind turbine generator can be calculated based on the following formula:

${\Delta{W_{PMSG}(t)}} = {{\int_{t}^{t + \frac{2\pi}{\omega}}{\Delta P_{e}d{\Delta\theta}_{pll}}} + {\int_{t}^{t + \frac{2\pi}{\omega}}\left( {{\Delta i_{d}d\Delta u_{q}} - {\Delta i_{q}d\Delta u_{d}}} \right)}}$

Where ΔP_(e) is the variation of direct-drive wind turbine generator output active power compared with the steady-state value; Δθ_(pll) is the variation of PLL angle; Δi_(d), Δi_(q), Δu_(q) and Δu_(d) are respectively the variation of d-axis and q-axis components of current and voltage at the terminal of direct-drive wind turbine generator, and ω is the oscillation frequency of dq axis.

Furthermore, the second processor assesses the stability level of system according to the value of η, including:

S1: discretizing the terminal energy:

${\Delta{W_{PMSG}(i)}} = {\sum\limits_{t = t_{i}}^{t_{i} + \frac{2\pi}{\omega}}{\Delta{{W_{PMSG}(t)}.}}}$

The value of terminal energy of direct-drive wind turbine generator at the time, t₁, . . . t_(n), i.e. ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) can be obtained; and n is the number of data;

S2: with discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)), calculating the energy negative gradient η, which includes the following steps in detail:

calculating the discretized average terminal energy of direct-drive wind turbine generator by the following formula:

${\Delta\;{\overset{\_}{W}}_{PMSG}} = {\frac{1}{n}{\sum\limits_{i = 1}^{n}{\Delta{{W_{PMSG}(i)}.}}}}$

Applying the discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) and its average value to the following formula, and obtaining the energy negative gradient based on it;

$\eta = {- {\nabla\left\lbrack \frac{{{\Delta\;{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}}}{\left( {\frac{1}{n}{\sum\limits_{t = 1}^{n}\left( {{\Delta{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}} \right)^{2}}} \right)^{\frac{1}{2}}} \right\rbrack}}$

Where ∇[⋅] is the symbol of gradient;

S3: assessing the stability level of system according to η: when η is positive, the system is stable and the bigger the value of η is, the higher the system stability level is; When η is zero, the system is critically stable; When η is negative, the system will go unstable, and as the absolute value of η increases, system stability level will drop.

Compared with the prior technology, the application can realize at least one of the following beneficial effects:

1. The calculation of energy negative gradient can quantify the stability level of the system, realize the early warning of system oscillation and ensure the safe and stable operation of the power grid.

2. Through the terminal energy calculation module, the calculation amount can be reduced, the occurrence and development of oscillation can be analyzed in real time, the calculation error can be reduced, and the accurate operation of the system can be guaranteed.

3. The parameter analysis module can give the influence of power grid and generator side parameters on the oscillation stability level, realize the mutual cooperation between the stability level evaluation and parameter adjustment, and adjust the parameters based on the adjustment rules, so as to make the system stability meet the requirements and improve the working efficiency of the system.

4. Through the simulation experiment, the stability evaluation and parameter analysis under different oscillation scenarios are realized, which improves the security and reliability of the system.

In the application, the above technical schemes can be combined with each other to realize more preferred combination schemes. Other features and advantages of the application will be described in a subsequent specification, and some of the advantages may become apparent from the specification or be understood by implementing the application. The object and other advantages of the application can be realized and obtained through the contents specially pointed out in the specification and the drawings.

BRIEF DESCRIPTION OF DRAWINGS

The attached figures are only for the purpose of illustrating specific embodiments, and are not considered to be a limitation of the application. In the whole figures, the same reference symbols represent the same parts.

FIG. 1 is the structure diagram of direct-drive wind turbine generator integrated system in an embodiment;

FIG. 2 is the structure diagram of a stability evaluation system of direct-drive wind turbine generator in another embodiment;

FIG. 3 is the simulation structure diagram of a stability evaluation system of direct-drive wind turbine generator in another embodiment;

FIG. 4 is time-domain simulation diagram of different wind turbines at stable operating points in another embodiment;

FIG. 5 is time-domain simulation diagram for different PLL parameters in another embodiment;

FIG. 6 is time-domain simulation diagram under different grid parameters in another embodiment;

FIG. 7 is the diagram of the energy negative gradient of direct-drive wind turbine generator at different stable operating points in another embodiment;

FIG. 8 is the change diagram of energy negative gradient components of direct-drive wind turbine generator at different stable operating points in another embodiment;

FIG. 9 is diagram of the energy negative gradient of direct-drive wind turbine generator for different PLL parameters in another embodiment;

FIG. 10 is the change diagram of energy negative gradient components of direct-drive wind turbine generator for different PLL parameters in another embodiment;

FIG. 11 is diagram of the energy negative gradient of direct-drive wind turbine generator for different grid lines parameters in another embodiment;

FIG. 12 is the change diagram of energy negative gradient components of direct-drive wind turbine generator for different grid lines parameters in another embodiment.

DESCRIPTION OF EMBODIMENTS

The preferred embodiments of the application will be described below in combination with the attached figures, where the attached figures form part of the application and, together with the embodiments of the application, are used to explain the principles of the application, not to define the scope of the application.

One embodiment of the application reveals a stability evaluation method of direct-drive wind turbine generator, including the following steps:

measuring the variation of voltage, current, active power and PLL (phase-looked-loop) angle at the terminal of direct-drive wind turbine generator in one oscillation period

$\frac{2\pi}{\omega};$

based on the measurements of the variation, calculating the terminal energy of direct-drive wind turbine generator with the established terminal energy model of direct-drive wind turbine generator;

according to the calculated terminal energy of direct-drive wind turbine generator, obtaining the energy negative gradient of direct-drive wind turbine generator, and assessing the stable level of the system according to the value of energy negative gradient; and

based on the influence pattern of the critical parameters of stable operation point of direct-drive wind turbine generator, PLL and transmission lines on the energy negative gradient, obtaining the adjustment law of the critical parameters that affect the stability of the system.

By calculating the energy negative gradient, the stable level of the system is quantified, the system oscillation warning is realized, and the safe and stable operation of the power grid is ensured. When the system stability level does not meet the requirements, parameters can be adjusted based on the adjustment law until the system stability meets the requirements.

Specifically, the structure of the system with direct-drive wind turbine generators integrated is shown in FIG. 1 . Establish the terminal energy model of direct-drive wind turbine generators with the critical parameters of PLL and transmission lines, which is:

${\Delta{W_{PMSG}(t)}} = {{\int_{t}^{t + \frac{2\pi}{\omega}}{\Delta\; P_{e}d\;\Delta\;\theta_{pll}}} + {\int_{t}^{t + \frac{2\pi}{\omega}}{\left( {{\Delta i_{d}d\Delta u_{q}} - {\Delta i_{q}d\Delta_{u_{d}}}} \right).}}}$

Where ΔP_(e) is the variation of direct-drive wind turbine generator output active power of direct-drive wind turbine generator compared with the stable value, Δθ_(pll) is the variation of PLL angle, Δi_(d), Δi_(q), Δu_(q) and Δu_(d) are respectively the variation of d-axis and q-axis components of current and voltage at the terminal of direct-drive wind turbine generator, and ω is the oscillation frequency of dq axis.

By monitoring the terminal energy of direct-drive wind turbine generator, it can analyze the occurrence and development of oscillation in real time and reduce the calculation amount and calculation error. The calculation efficiency is high and it ensures the accurate operation of the system.

Preferably, the calculation process of energy negative gradient includes:

step 1: discretizing the terminal energy:

${\Delta{W_{PMSG}(i)}} = {\sum\limits_{t = t_{i}}^{t_{i} + \frac{2\;\pi}{\omega}}{\Delta\;{{W_{PMSG}(t)}.}}}$

The value of terminal energy of direct-drive wind turbine generator at the time t₁, . . . t_(n), i.e. ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) can be calculated; n is the number of data.

Step 2: with discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) to calculate the energy negative gradient η, which includes the following steps in detail.

Calculating the average value of the discretized terminal energy of direct-drive wind turbine generator by the following formula:

${\Delta\;{\overset{\_}{W}}_{PMSG}} = {\frac{1}{n}{\sum\limits_{i = 1}^{n}{\Delta{{W_{PMSG}(i)}.}}}}$

With the discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) and the average value, the energy negative gradient can be obtained based on the following formula.

$\eta = {- {\nabla\left\lbrack \frac{{{\Delta\;{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}}}{\left( {\frac{1}{n}{\sum\limits_{t = 1}^{n}\left( {{\Delta{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}} \right)^{2}}} \right)^{\frac{1}{2}}} \right\rbrack}}$

Where ∇[⋅] is the symbol of gradient.

Step 3: assessing the stability according to η: when η>0, the system is stable, and the bigger the value of η is, the higher the system stability level is. when η=0, the system is critically stable. when η<0, the system will go unstable, and as the absolute value of η increases, the system stability level will drop.

η being positive or negative respectively represents the ascending or descending trend of dynamic energy of direct-drive wind turbine generator, i.e. the system being stable or unstable. The absolute value of η represents the ascending or descending speed of dynamic energy. When η is positive, the bigger the absolute value of η is, the higher the descending amount of dynamic energy per unit time is, and the faster system oscillation converges, the higher system stability level is. Conversely, the smaller the absolute value of η is, the more slowly system oscillation converges, and the lower system stability level is. When η is negative, the system will go unstable, and as the absolute value of η increases, the system stability level will drop and system oscillation will deteriorate. Therefore, the stable level of sub/super synchronous oscillation of the system can be assessed quantitatively according to the energy negative gradient η.

Through the energy negative gradient of direct-drive wind turbine generator, the stability of the system can be judged more quickly and effectively, the real-time monitoring of dynamic energy can be realized, and the working efficiency of the system can be improved.

Preferably, the adjustment law of the critical parameters affecting the stability of the system include:

step S1: applying the PLL dynamic function to the energy negative gradient expression, the detailed expression with critical parameters of stable operating points and the PLL can be achieved:

$\eta_{2} = {{- \left( {{I_{{sd}\; 0}\omega\frac{\cos\;\varphi_{c}}{2}} + {I_{{sq}\; 0}\omega\frac{\sin\;\varphi_{c}}{2}}} \right)}{\quad{{{\quad\left\lbrack {{{A_{-}(t)}\Omega_{u -}\alpha_{u}e^{\alpha_{u}t}} + {\frac{{dA}_{-}(t)}{dt}\Omega_{u -}e^{\alpha_{u}t}}} \right\rbrack }_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}} - {U_{{sd}\; 0}{\quad{\omega{\frac{\cos\;\left( {\varphi_{c} - \varphi_{RL}} \right)}{2}\left\lbrack {{{A_{-}(t)}\alpha_{i}\Omega_{i -}\ e^{\alpha_{i}t}} + {\frac{{dA}_{-}(t)}{dt}\Omega_{i -}\ e^{\alpha_{i}t}}} \right\rbrack}}}_{t_{0}}^{t_{0} + \frac{2\;\pi}{\omega}}} + {\left( {{I_{{sd}\; 0}\omega\frac{\cos\;\varphi_{c}}{2}} + {I_{{sq}\; 0}\omega\frac{\sin\;\varphi_{c}}{2}}} \right)\left\lbrack {{{A_{+}(t)}\Omega_{u +}\alpha_{u}e^{\alpha_{u}t}} + {\frac{d{A_{+}(t)}}{dt}\Omega_{u +}e^{\alpha_{u}t}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}{{+ U_{{sd}\; 0}}\omega{\frac{\cos\left( {\varphi_{c} - {\varphi\phi}_{RL}} \right)}{2}\left\lbrack {{{A_{+}(t)}\alpha_{i}\Omega_{i +}e^{\alpha_{i}t}} + {\frac{d{A_{+}(t)}}{dt}\Omega_{i +}e^{\alpha_{i}t}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}{{- 2}\;\omega\; K_{pll}^{2}\Omega_{u -}\Omega_{u +}e^{2\alpha_{u}t}\alpha_{u}\sin\;{\varphi_{c}\left\lbrack {{I_{sd0}{\sin\left( {\varphi_{u +} + \varphi_{u -}} \right)}} - {I_{sq0}{\cos\left( {\varphi_{u +} + \varphi_{u -}} \right)}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}}}}$

Where A⁻, A₊ respectively represent the amplitude of PLL angle variation Δθ_(pll) in sub-synchronous frequency and super-synchronous frequency. I_(sd0), I_(sq0), U_(sd0), U_(sq0) are the original value of dq axes components of current and voltage of the direct-drive wind turbine generator terminal. α_(u) and α_(i) are the attenuation coefficients of oscillation components of voltage and current. Ω_(u−), Ω_(u+), Ω_(i−) and Ω_(i+) are respectively the amplitude coefficients of sub/super synchronous components of voltage and current. φ_(u−) and φ_(u+) are the initial phase angles of sub-synchronous and super-synchronous voltages. φ_(i−) and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents. Ω_(u+)=Ω_(i+)X₊, Ω_(u−)=Ω_(i−)X⁻, φ_(u−)=φ_(i−)+φ_(RL) and φ_(u+)=φ_(i+)+φ_(RL), where φ_(RL) is the phase angle of transmission lines. k_(pθ) and k_(iθ) are the proportion gain and integration gain of PLL.

${\varphi_{c} = {\arctan\frac{\omega k_{p\theta}}{k_{i\;\theta}\;}}},{K_{pll} = {{\frac{{k_{p\theta}s} + k_{i\;\theta}}{s^{2} - {U_{sq0}\left( {{k_{p\theta}s} + k_{i\;\theta}} \right)}}}_{s = {j\;\omega}}.}}$

Applying the dynamic function of transmission lines to expression of the energy negative gradient, the one with transmission lines parameters can be expressed as:

${{{\eta_{1} = {{- 2}\Omega_{i +}\Omega_{i -}e^{2\alpha_{i}t}\alpha_{i}X_{L}\omega\;{\sin\left( {\varphi_{i +} + \varphi_{i -}} \right)}}}}_{t_{0}}^{t_{0} + \frac{2\;\pi}{\omega}} + {2\left\lbrack {{\Omega_{i -}^{2}e^{2\alpha_{i}t}\alpha_{i}\omega\; R} - {\Omega_{i +}\Omega_{i -}e^{2\alpha_{i}t}\alpha_{i}R\;{{\omega cos}\left( {\varphi_{i +} + \varphi_{i -}} \right)}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}.$

Where Ω_(i−) and Ω_(i+) are respectively the amplitude coefficients of sub-synchronous component and super-synchronous component of current. α_(i) is the attenuation coefficients of oscillation components of current. R and X_(L) are line resistance and the equivalent reactance of the transmission network. φ_(i−) and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents.

Step S2: based on the constructed energy negative gradient with the critical parameters of stable operating points and PLL and the one with the critical parameters of the transmission lines, the influence pattern of these parameters on the energy negative gradient can be analyzed.

The analysis process of the influence pattern includes:

The first two items of the above energy negative gradient model with critical parameters of transmission lines are relevant to only line reactance X_(L). Since

${{e^{2\;\alpha_{i}t}\alpha_{i}}}_{t_{0}}^{t_{0} + \frac{2\;\pi}{\omega}}$ is constantly above 0, this component of line energy negative gradient is constantly below 0. According to the energy negative gradient model with grid parameters, the bigger X_(L) is, the weaker power grid strength is, the smaller line energy negative gradient is, and the lower the system stability level is.

The last two items of the above energy negative gradient model with critical parameters of transmission lines are relevant to only line parameter R, and the value of this component is affected by the amplitude coefficients of sub/super synchronous currents. When sub-synchronous component of the system is dominant, i.e. when Ω_(i−)>Ω_(i+), the component of line energy negative gradient corresponding to the resistance R is constantly positive, and the bigger the resistance R is, the bigger line energy negative gradient is. When super-synchronous component is dominant, the energy negative gradient generated by the resistance R can be determined by the value (that is, plus or minus sign of the value) of Ω_(i−)−Ω_(i+) cos(φ_(i+)+φ_(i−)). In order to obtain the values of Ω_(i−) and Ω_(i+) cos(φ_(i+)+φ_(i−)), the expression of cos(φ_(i+)+φ_(i−)) is further analyzed.

${\varphi_{i +} + \varphi_{i -}} \approx {{\Delta\theta}_{pll} + \varphi_{RL} + {\frac{\pi}{2}.}}$

Where φ_(i−) and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents. Δθ_(pll) is the variation of PLL angle. φ_(RL) is the impedance angle of transmission lines.

The value of cos(φ_(i+)+φ_(i−)) is mainly determined by Δθ_(pll) and φ_(RL). On the one hand, Δθ_(pll) is mainly affected by PLL parameters. Currently, PLL usually applies fixed voltage phase control, thus K_(pθ)<<K_(iθ) and the value of Δθ_(pll) is relatively small during sub/super synchronous oscillation. On the other hand, in order to reduce the transmission loss on transmission lines, the value of transmission lines resistance R is usually small, especially in long-distance transmission cases R<<X_(L). Therefore, Δθ_(pll)+φ_(RL)≈0 and the value of cos(φ_(i+)+φ_(i−)) is approximately 0 and constantly smaller than

$\frac{\Omega_{i -}}{\Omega_{i +}}.$ Thus, according to the energy negative gradient model with grid parameters, the energy negative gradient of line generated by R is positive, and as R increases, the energy negative gradient also increases, and the stable level of the system improves.

Based on the above analysis, in two scenarios where sub-synchronous and super-synchronous is dominate, the energy negative gradient generated by line resistance is constantly positive and increases as R increases.

The first two items of the energy negative gradient model with PLL are energy negative gradient item being relevant to only sub-synchronous components and are denoted as η⁻; the following two items are energy negative gradient item being relevant to only super-synchronous components and are denoted as η₊; the last item is affected by both sub-synchronous and super-synchronous components, which represents the coupling effect between sub-synchronous components and super-synchronous components generated by PLL dynamic characteristics and is denoted as η_(−∩+). With the above three parts synthesized, the energy negative gradient model with PLL can be written as: η₂=−η⁻+η₊+η_(−∩+)

When sub-synchronous oscillation is dominant in the system, η₊−η⁻<0. Since direct-drive wind turbine generator operates normally under condition of the unity power factor i.e. I_(sd0)>>I_(sq0), η_(−∩+)<0, i.e. the coupled energy negative gradient generated by PLL is constantly negative. In this case, according to the energy negative gradient model with PLL, η₂ is constantly below 0, and the energy negative gradient generated by PLL is constantly negative. If PLL parameters K_(p) and K_(i) are increased, η₂ will decrease and the stable level of the system will drop.

When super-synchronous oscillation is dominant in the system, η₊−η⁻>0 and the coupled energy negative gradient generated by PLL η_(−∩+) is constantly negative. According to the energy negative gradient model with PLL, the value (that is, plus or minus sign of the value) of η₂ is determined by the values of η_(−∩+) and η₊−η⁻. When η_(−∩+)>η⁻−η₊, the energy negative gradient generated by PLL is positive; otherwise, the energy negative gradient generated by PLL is negative. PLL parameters affect the energy negative gradient by changing K_(pll) in A⁻ and A₊. If PLL parameters K_(p) and K_(i) are increased, K_(pll) will increase, η₊−η⁻ will also increase, but η_(−∩+) will decrease. Furthermore, according to the energy negative gradient model with PLL, η_(−∩+) is inversely proportional to K_(pll) ², and as PLL control parameters vary, the decreasing amplitude of η_(−∩+) is larger than the increasing amplitude of η₊−η⁻, thus η₂ as a whole exhibits decreasing trend. Based on the two scenarios of sub-synchronous oscillation or super-synchronous oscillation being dominant, increasing PLL parameters will lower the stable level of power system integrated with direct-drive wind turbine generator.

It can be seen from the energy negative gradient model with PLL that, the energy negative gradient of PLL is also affected by current parameters at the equilibrium point of wind turbine, i.e. I_(sd0) and I_(sq0), which are determined by the reference values of grid-side converter of direct-drive wind turbine generator.

When sub-synchronous component is dominant, η₂<0, and according to the energy negative gradient model of PLL, the energy negative gradient of PLL will decrease as I_(sd0) increases and increase as I_(sq0) increases.

When super-synchronous component is dominant, due to the limits of wind turbine LVRT (Low voltage ride through) and HVRT (High voltage ride through), the PLL bandwidth of direct-drive wind turbine generator cannot be set too low, thus the value of K_(pll) in super-synchronous frequency band is relatively large. In this case, according to the model of the energy negative gradient of PLL, η_(−∩+)<η⁻−η₊, i.e. η₂<0. If I_(sd0) is increased, the energy negative gradient of PLL will further decrease, aggravating the stable level of the system. On the other hand, if I_(sq0) is increased, η₊−η⁻ will also increase, the absolute value of η_(−∩+) will decrease, η₂ will increase, and the stable level of the system will be improved. Therefore, based on the two cases of sub-synchronous oscillation or super-synchronous oscillation being dominant respectively, moderately decreasing I_(sd0) and increasing I_(sq0) will improve the stable level of the system.

Based on the above analysis, the influence patterns are included: Reducing Isd0 or increasing Isq0 help improve the stable level of the system. The bigger R is and the smaller XL is, the more stable the system is. Reducing PLL parameters K_(p) and K_(i), the system stability can be enhanced.

Through parameter analysis, the influence of power grid and generator side parameters on the level of oscillation stability is given. The coordination of stable level evaluation and parameter adjustment is realized, and the working efficiency of the system is improved.

Another embodiment of the application is shown in FIG. 2 , which provides a stability evaluation system of direct-drive wind turbine generator, including a data acquisition module, a terminal energy calculation module, a stability evaluation module, a critical parameter analysis module and a result output module. Wherein the data acquisition is implemented by the PMU, the terminal energy calculation is implemented by a first processor, the stability evaluation is implemented by a second processor, critical parameter analysis is implemented by a third processor and the result output is implemented by the result output terminal. Each of the first processor, the second processor and the third processor is independent processor, or the three processors are integrated in a single processor. The data acquisition module is used to collect the variation of voltage, current, active power and PLL angles at the terminal of direct-drive wind turbine generator. The terminal energy calculation module is used to calculate the terminal energy of direct-drive wind turbine generator. The stability evaluation module is used to assess the stable level of direct-drive wind turbine generator according to the value of energy negative gradient. The critical parameter analysis module is used to analyze the influence patterns of stable operation point of direct-drive wind turbine generator, PLL and transmission lines parameters on the energy negative gradient and propose the adjustment law of critical parameters to improve the stability of the system. The result output module is used to output the stability level of direct-drive wind turbine generator and adjustment law of the critical parameters affecting system stability.

By measuring the variation of terminal voltage, current, active power and locked phase angle at the terminal of direct-drive wind turbine generator, the system can effectively monitor the stable level of sub/super synchronous oscillations in the grid connected system of direct-drive wind turbine generator, and reveal the influence pattern of system parameters on the stable level. Based on the influence law, the parameters can be adjusted to improve the stability and safety of the system.

Optimally, the data acquisition module collects the variation of voltage, current, active power and PLL angles by the PMU (Phasor Measurement Unit) equipped at the terminal of direct-drive wind turbine generator and PLL. Through the data acquisition module, the comprehensiveness and accuracy of data acquisition are realized.

Optimally, the terminal energy calculation module receives the data transferred by data acquisition module. And the terminal energy of direct-drive wind turbine generator can be calculated:

${\Delta{W_{PMSG}(t)}} = {{\int_{t}^{t + \frac{2\pi}{\omega}}{\Delta\; P_{e}d\;\Delta\;\theta_{pll}}} + {\int_{t}^{t + \frac{2\pi}{\omega}}\left( {{\Delta i_{d}d\Delta u_{q}} - {\Delta i_{q}d\Delta_{u_{d}}}} \right)}}$

Where ΔP_(e) is the variation of direct-drive wind turbine generator output active power compared with the stable value. Δθ_(pll) is the variation of PLL angle. Δi_(d), Δi_(q), Δu_(q) and Δu_(d) are respectively the variation of d-axis and q-axis components of current and voltage at the terminal of direct-drive wind turbine generator. ω is the oscillation frequency of dq axis.

Through the terminal energy calculation module, the calculation amount is reduced, the occurrence and development of oscillation can be analyzed in real time, the calculation error is reduced, and the accurate operation of the system is guaranteed.

Optimally, the stability evaluation module can assess the stable level of the system according to η, including:

step S1: discretizing the terminal energy:

${\Delta{W_{PMSG}(i)}} = {\sum\limits_{t = t_{i}}^{t_{i} + \frac{2\;\pi}{\omega}}{\Delta\;{{W_{PMSG}(t)}.}}}$

The value of terminal energy of direct-drive wind turbine generator at the time, t₁, . . . t_(n), i.e. ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) can be calculated. n is the number of data.

Step S2: with discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)), calculating the energy negative gradient η, which includes the following steps.

Calculating the average of the discretized terminal energy of direct-drive wind turbine generator by the following formula:

${\Delta\;{\overset{\_}{W}}_{PMSG}} = {\frac{1}{n}{\sum\limits_{i = 1}^{n}{\Delta{{W_{PMSG}(i)}.}}}}$

With the discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) and its average value, the energy negative gradient can be obtained by the following formula.

$\eta = {- {\nabla\left\lbrack \frac{{{\Delta\;{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}}}{\left( {\frac{1}{n}{\sum\limits_{t = 1}^{n}\left( {{\Delta{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}} \right)^{2}}} \right)^{\frac{1}{2}}} \right\rbrack}}$

Where ∇[□] is the symbol of gradient.

Step S3: the stable level of system can be assessed according to η. When η is positive, the system is stable and the bigger the value of η is, the higher system stability level is. When η=0, the system is critically stable. When η is negative, the system will go unstable, and as the absolute value of η increases, system stability level will drop.

The stability evaluation module can judge the stability of the system quickly and accurately, improve the working efficiency of the system, and ensure the stable operation of the system.

To verify the feasibility of the above proposed method, the simulation system integrated with direct-drive wind turbine generators has been constructed. With direct-drive wind turbine generators operating under the condition of different PLL parameters, transmission lines parameters, the stable level of system can be assessed, corresponding parameters can be adjusted so that the system is stable, and the stability assessment results are verified by simulation experiments.

Specifically, the simulation structure of the system is shown in FIG. 3 . The model of power system integrated with direct-drive wind turbine generators (it is PMSG (Permanent Magnetic Synchronous Generator) in FIG. 3 ) is constructed in RT-LAB. Parameter settings of power system integrated with direct-drive wind turbine generator are shown in Table 1:

TABLE 1 Parameter settings of power system integrated with direct-drive wind farm Parameter Symbol Value PMSG current inner loop proportion gain k_(p1)/p.u.   0.201 PMSG current inner loop integration gain k_(i1)/p.u. 40.2 PMSG PLL proportion gain k_(pθ)/p.u.  1  PMSG PLL integration gain k_(iθ)/p.u.  40  Amplitude of sub synchronous Ω_(u−)/p.u.   0.95 component of voltage at PCC Amplitude of sub synchronous Ω_(i−)/p.u.    1.0834 component of current at PCC Amplitude of super synchronous Ω_(u+)/p.u.    0.969 component of voltage at PCC Amplitude of super synchronous Ω_(i+)/p.u.    1.107 component of current at PCC Phase of sub synchronous component φ_(U−)/deg    −0.0125 of voltage at PCC Phase of sub synchronous component φ_(I−)/deg    19.2779 of current at PCC Filter inducance L/H  2e−3

Considering different operation points, PLL and transmission lines parameters, the operating results of stability assessment module are shown as follows.

Three groups of parameters, i.e. Isq0=1, Isd0=1.3, Isq0=1.3, Isd0=1.3; Isq0=1.3, Isd0=1 are selected for time-domain simulation tests, and the corresponding calculated values of energy negative gradient are shown in Table 2:

TABLE 2 Parameter settings of time-domain simulation and calculated results of energy negative gradient energy negative Oscillation Parameters gradient characteristic Isq0 = 1, Isd0 = 1.3 −1.02 Diverging to instability Isq0 = 1.3, Isd0 = 1.3 −0.73 Diverging slowly Isq0 = 1.3, Isd0 = 1   0.0503 Converging to stability

Time-domain simulation of direct-drive wind turbine generator at different stable operating points are shown in FIG. 4 . It can be seen from FIG. 4(a) that, when Isq0 increases from 1 to 1.3, the amplitude of system oscillation obviously decreases, and the stable level of the system improves. However, in the case of Isq0=1.3, Isd0=1.3, system oscillation still exhibits diverging trend, and the stable level of the system is still negative, which is consistent with the calculated value of energy negative gradient shown in Table 2. According to FIG. 4(b), when Isd0 increases from 1 to 1.3, system oscillation turns from converging to diverging, and system stability level drops acutely. Based on the analysis of calculated values in Table 2, the value of energy negative gradient decreases from 0.0503 to −0.73, the system turns from converging rapidly to diverging slowly, which is consistent with the stability level analysis results from time domain simulation.

Four groups of PLL parameter settings are selected for time-domain simulation tests. The settings of parameters and the corresponding calculated values of energy negative gradient are shown in Table 3.

TABLE 3 Parameter settings of time-domain simulation and calculated results of energy negative gradient energy negative Oscillation Parameters gradient characteristic Kp = 1.85, Ki = 100 −0.9  Diverging acutely Kp = 1.15, Ki = 100  −0.0723 Diverging slowly Kp = 1.2, Ki = 125 −0.92 Diverging acutely Kp = 1.2, Ki = 55 −0.04 Diverging slowly

Time-domain simulation for different PLL parameters are shown in FIG. 5 . FIG. 5(a) shows the time-domain simulation curves of d-axis current when Ki=100, Kp=1.85 and Ki=100, Kp=1.15. It can be seen that, as Kp decreases from 1.85 to 1.15, the amplitude of system oscillation obviously drops, the diverging trend slows down, and system stability level improves. In Table 3, the calculated values of energy negative gradient corresponding to parameter settings of Kp=1.85, Ki=100, Kp=1.15 are respectively−0.9 and −0.0723, with stability analysis results being diverging acutely and diverging slowly, which are consistent with the variation trends of time-domain simulation curves. When Kp=1.2, Ki=125 and Kp=1.2, Ki=55, time-domain simulation curves of d-axis current are shown in FIG. 5(b). It can be seen that, as Ki decreases, diverging of system oscillation obviously slows down, and system stability level improves. In Table 3, the calculated values of energy negative gradient corresponding to parameter settings of Kp=1.2, Ki=125 and Kp=1.2, Ki=55 are respectively−0.92 and −0.04, which are also consistent with the diverging trends of time-domain simulation curves.

Four groups of grid parameter settings are selected for time-domain simulation tests. The settings of parameters and the corresponding calculated values of energy negative gradient are shown in Table 4.

TABLE 4 Parameter settings of time-domain simulation and calculated results of energy negative gradient energy negative Oscillation Parameters gradient characteristic R = 0.018, XL = 1.13 −0.2526 Diverging to instability R = 0.15, XL = 1.13  0.04271 Converging to stability R = 0.15, XL = 1  0.02364 Converging to stability R = 0.15, XL = 2.5 −0.1764 Diverging to instability

Time-domain simulation diagrams under different grid parameters are shown in FIG. 6 . It can be seen from FIG. 6(a) that, when line resistance R increases from 0.018 p.u. to 0.15 p.u., system oscillation turns from diverging acutely to converging obviously, and system stability level improves. the calculated values of energy negative gradient corresponding to the first two groups of parameters in Table 4, are respectively−0.2526 and 0.04271, which are consistent with the variation trends of time-domain simulation curves. According to FIG. 6(b), when grid equivalent reactance increases from 1 p.u. to 2.5 p.u., system oscillation turns from converging to diverging, and system stability level drops acutely. the calculated values of energy negative gradient corresponding to the last two groups of parameters in Table 4 are respectively 0.02364 and −0.1764, which are consistent with the variation trends of time-domain simulation results.

Critical parameters analysis module analyzes the influence of different parameters on system stability. The operating results are shown as follows.

When Isd0 and Isq0 increase from 0.16 to 1.67 in the step of 0.08, the variation trend of the energy negative gradient is depicted, as shown in FIG. 7 . It can be seen that, the energy negative gradient of system is mainly affected by Isd0. As Isd0 increases, the energy negative gradient drops significantly. The effect of Isq0 on the energy negative gradient of system is relatively small. As Isq0 increases, the energy negative gradient of system slightly increases, and system stability level improves.

The variation trends of components of energy negative gradient corresponding to different stable operating points of direct-drive wind turbine generator are shown in FIG. 8 . It can be seen that, when Isq0=1.08 and Isd0 increases, and the energy negative gradient of lines scarcely varies, while the energy negative gradient of PLL decreases, gradually offsetting the positive dissipation effect of transmission line and causing system stability level to drop. When Isd0=1.08 and Isq0 increases, the energy negative gradient of PLL increases, and system stability level gradually improves.

With the other system operation parameters remaining unchanged, PLL parameter Kp increases from 0.1 to 2 in the step of 0.05, and parameter Ki increase from 10 to 200 in the step of 5, and the energy negative gradient corresponding to different parameters can be calculated, the results is shown in FIG. 9 . It can be seen that, when Kp<1 and Ki<115, the energy negative gradient of system is positive, i.e. the system is stable. As Kp and Ki gradually increase, the energy negative gradient of the system decreases.

FIG. 10 shows the variation of two components of energy negative gradient of direct-drive wind turbine generator corresponding to different values of PLL parameters. When Ki=100 and the value of Kp varies, the variation curves of energy negative gradient of lines and that of PLL are shown in FIG. 10(a). As shown in figure, in the one hand, the energy negative gradient of line parameters is positive and scarcely affected by PLL parameters; in the other hand, the energy negative gradient of PLL is negative and gradually decreases as Kp increases, the positive dissipation effect of transmission lines can be offset and the overall stable level of system can be lowered. When Kp=1, energy negative gradient varies with distinct variation of Ki, the variation curves of energy negative gradient of lines and that of PLL are shown in FIG. 10(b). It can be seen that, the energy negative gradient of lines is positive and gradually decreases as Ki increases. This is because the proportion of super-synchronous component is raised when PLL parameter increases, the energy negative gradient of lines is caused to drop. On the other hand, the energy negative gradient of PLL is negative and decreases as Ki increases. The decreasing amplitude of energy negative gradient of PLL is much larger than that of energy negative gradient of lines, thus the energy negative gradient of system overall exhibits descending trend, and stability level of the system drops as Ki increases.

FIG. 11 shows the effect of different transmission lines parameters on the energy negative gradient. where line resistance R increases from 0.01 p.u. to 0.2 p.u. in the step of 0.005 p.u., and grid equivalent reactance XL increases from 0.5 p.u. to 3.5 p.u. in the step of 0.05 p.u. It can be seen that, as line resistance increases, the energy negative gradient of system increasingly sharply. The effect of grid equivalent reactance on the energy negative gradient is relatively smaller than the effect of line resistance on it. As grid equivalent reactance increases, i.e. as the strength of grid decreases, the energy negative gradient of system also decreases, and system stability level drops.

The variation curves of components of energy negative gradient corresponding to different values of grid parameters are shown in FIG. 12 . FIG. 12(a) shows the variation trends of energy negative gradient of lines and that of PLL when grid equivalent reactance remains 1.06 p.u. and line resistance increases. The energy negative gradient of PLL is relevant to line impedance angle. Since line resistance varies between 0.01 p.u. and 0.2 p.u., line impedance angle only varies very small in a small range between 79.3° and 88.46°. Thus, the energy negative gradient of PLL is scarcely affected by line resistance, while the energy negative gradient of line is almost linearly correlated with line resistance and increases significantly as line resistance increases. FIG. 12(b) shows the variation curves of components of energy negative gradient when line resistance remains 0.04p.u. and the strength of grid decreases. It can be seen that, the effect of grid strength on the energy negative gradient of line is much smaller than the effect of line resistance on it. However, according to the enlarged display, the energy negative gradient of line decreases as the strength of grid decreases.

A non-transitory machine-readable storage medium comprising instructions that when executed cause a processor of a computing device to: collect variation of voltage, current, active power and PLL angle at the terminal of direct-drive wind turbine generator; calculate the terminal energy of direct-drive wind turbine generator based on the variation; assess the stability level of direct-drive wind turbine generator according to the value of energy negative gradient; analyze the influence patterns of stable operation points, PLL and transmission lines parameters on the energy negative gradient and propose the adjustment law of critical parameters to improve the stability of the system; and output the stability level of the system and adjustment law of the critical parameters affecting system stability.

Through the simulation experiment, the stability evaluation and parameter analysis under different oscillation scenarios are realized, which improves the security and reliability of the system.

It can be understood by those skilled in the art that all or part of the flow of the method of the above embodiment can be realized by instructing related hardware through a computer program, which can be stored in a computer-readable storage medium. Wherein, the computer-readable storage medium is a disk, an optical disc, a read-only storage memory, a random storage memory, and the like.

The above is only a preferred embodiment of the application, but the scope of protection of the application is not limited to this. Any change or replacement that can be easily thought of by any person familiar with the technical field within the technical scope of the application shall be covered in the scope of protection of the application

The foregoing descriptions of specific exemplary embodiments of the present application have been presented for purposes of illustration and description. They are not intended to be exhaustive or to limit the application to the precise forms disclosed, and obviously many modifications and variations are possible in light of the above teachings. The exemplary embodiments were chosen and described in order to explain certain principles of the application and their practical application, to thereby enable others skilled in the art to make and utilize various exemplary embodiments of the present application, as well as various alternatives and modifications thereof. It is intended that the scope of the application be defined by the Claims appended hereto and their equivalents. 

What is claimed is:
 1. A stability evaluation method of direct-drive wind turbine generator, comprising: measuring, with a Phasor Measurement Unit (PMU), variations of voltage, current, active power and PLL (Phase Locked Loop) angle at a terminal of direct-drive wind turbine generator in one oscillation period $\frac{2\pi}{\omega};$ based on the measured variations, calculating terminal energy of direct-drive wind turbine generator with a terminal energy model of direct-drive wind turbine generator, wherein the terminal energy model of direct-drive wind turbine generator comprises: ${\Delta{W_{PMSG}(t)}} = {{\int_{t}^{t + \frac{2\pi}{\omega}}{\Delta\; P_{e}d\;\Delta\;\theta_{pll}}} + {\int_{t}^{t + \frac{2\pi}{\omega}}\left( {{\Delta i_{d}d\Delta u_{q}} - {\Delta i_{q}d\Delta u_{d}}} \right)}}$ where ΔP_(e) is the variation of direct-drive wind turbine generator output active power compared with the steady-state value, Δθ_(pll) is the variation of PLL angle, Δi_(d), Δi_(q), Δu_(q) and Δu_(d) are respectively the variation of d-axis and q-axis components of current and voltage at the terminal of direct-drive wind turbine generator, ω is the oscillation frequency of dq axis; obtaining energy negative gradient of direct-drive wind turbine generator according to discretized and average values of the terminal energy of direct-drive wind turbine generator, and assessing system stability level according to a value of the energy negative gradient, wherein the larger the value of the energy negative gradient is, the higher the system stability level is; and obtaining an adjustment law of critical parameters affecting the stability of the system based on the energy negative gradient with stability operation points and critical parameters PLL and transmission lines, which comprises the following steps S1-S2: step S1: applying the PLL dynamic function to the energy negative gradient expression, achieving the energy negative gradient with stability operation points and PLL critical parameters: $\left. {\eta_{2} = {{- \left( {{I_{sd0}\omega\frac{\cos\varphi_{c}}{2}} + {I_{sq0}\omega}} \right)}\frac{\sin\varphi_{c}}{2}}} \right)$ $\left\lbrack {{{A\_}(t)\Omega_{u -}\alpha_{u}e^{\alpha_{u}t}} + {\frac{{dA\_}(t)}{dt}\Omega_{u -}e^{\alpha_{u}t}}} \right\rbrack ❘_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}} -$ ${\left. {U_{sd0}\omega} \right){\frac{\cos\left( {\varphi_{c} - \varphi_{RL}} \right)}{2}\left\lbrack {{{A\_}(t)\alpha_{i}\Omega_{i -}e^{\alpha_{i}t}} + {\frac{{dA\_}(t)}{dt}\Omega_{i -}e^{\alpha_{i}t}}} \right\rbrack}❘_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}} +$ $\left. {{\left( {I_{sd0}\omega} \right)\frac{\cos\varphi_{c}}{2}} + {I_{sq0}\omega\frac{\sin\varphi_{c}}{2}}} \right)\left\lbrack {{{A_{+}(t)}\Omega_{u +}\alpha_{u}e^{\alpha_{u}t}} +} \right.$ $\left. {\frac{{dA}_{+}(t)}{dt}\Omega_{u +}e^{\alpha_{u}t}} \right\rbrack|_{t_{0}}^{{t_{0} + \frac{2\pi}{\omega} + {U_{sd0}\omega{\frac{\cos({\varphi_{c} - \varphi_{RL}})}{2}\lbrack{{{A_{+}(t)}\alpha_{i}\Omega_{i +}e^{\alpha_{i}t}} + {\frac{{dA}_{+}(t)}{dt}\Omega_{i +}e^{\alpha_{i}t}}}\rbrack}}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}} -$ ${2\omega K_{pll}^{2}\Omega_{u -}\Omega_{u +}e^{2\alpha_{u}t}\alpha_{u}\sin{\varphi_{c}\left\lbrack {{I_{sd0}{\sin\left( {\varphi_{u +} + \varphi_{u}} \right)}} - {I_{sq0}{\cos\left( {\varphi_{u +} + \varphi_{u -}} \right)}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}$ Where A⁻, A₊ respectively represent the amplitude of PLL angle variation Δθ_(pll) in sub-synchronous frequency and super-synchronous frequency; I_(sd0), I_(sq0), U_(sd0), U_(sq0) are the original value of dq axes components of current and voltage of the terminal of direct-drive wind turbine generator; α_(u) and α_(i) are the attenuation coefficients of oscillation components of voltage and current; Ω_(u−), Ω_(u+), Ω_(i−) and Ω_(i+) are respectively the amplitude coefficients of sub/super synchronous components of voltage and current; φ_(u−) and φ_(u+) are the initial phase angles of sub-synchronous and super-synchronous voltage; φ_(i−) and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents; Ω_(u+)=Ω_(i+)X₊, Ω_(u−)=Ω_(i−)X⁻, φ_(u−)=φ_(i−)+φ_(RL) and φ_(u+)=φ_(i+)+φ_(RL), where φ_(RL) is the phase angle of transmission lines; k_(pθ) and k_(iθ) are the proportion gain and integration gain of PLL; ${\varphi_{c} = {\arctan\frac{\omega k_{p\theta}}{k_{i\theta}}}},{{K_{pll} = {❘\frac{{k_{p\theta}s} + k_{i\theta}}{s^{2} - {U_{sq0}\left( {{k_{p\theta}s} + k_{i\theta}} \right)}}❘}_{s = {j\omega}}};}$ applying the dynamic function of system transmission lines to the energy negative gradient expression, obtaining the energy negative gradient of critical parameters of transmission lines: ${{\eta_{1} = {{- 2}\Omega_{i +}\Omega_{i -}e^{{2\alpha},t}\alpha_{i}X_{L}{{\omega sin}\left( {\varphi_{i +} + \varphi_{i -}} \right)}}}❘}_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}} + {{2\left\lbrack {{\Omega_{i -}^{2}e^{{2\alpha},t}\alpha_{i}\omega R} - {\Omega_{i +}\Omega_{i -}\ e^{{2\alpha},t}\alpha_{i}R{{\omega cos}\left( {\varphi_{i +} + \varphi_{i -}} \right)}}} \right\rbrack}❘_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}}$ where Ω_(i−) and Ω_(i+) are the amplitude coefficients of current sub-synchronous component and super-synchronous component respectively; α_(i) is the attenuation coefficients of oscillation components of current; R and X_(L) are line resistance and the equivalent reactance of the transmission lines; φ¹⁻ and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents; and step S2: based on the constructed energy negative gradient with the critical parameters of stability operation points of direct-drive wind turbine generator, PLL and the transmission lines, obtaining influence pattern of these parameters on the energy negative gradient; wherein the adjustment law comprises: the stability level of system is improved by reducing I_(sd0) or increasing I_(sq0); the more stable the system is the larger R is and the smaller X_(L) is; the system stability is enhanced by reducing PLL parameters K_(p) and K_(i).
 2. The stability evaluation method of direct-drive wind turbine generator according to claim 1, the calculation process of energy negative gradient comprises: step 1: discretizing the terminal energy: ${\Delta{W_{PMSG}(i)}} = {\sum\limits_{t = t_{i}}^{t + \frac{2\pi}{\omega}}{\Delta{W_{PMSG}(t)}}}$ the value of terminal energy of direct-drive wind turbine generator at the time t₁, . . . t_(n), i.e. ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) can be obtained; n is the number of data; and step 2: with discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)), calculating the energy negative gradient η.
 3. The stability evaluation method of direct-drive wind turbine generator according to claim 2, the step 2 comprises: calculating the discretized average value of the discretized terminal energy of direct-drive wind turbine generator by the following formula; ${\Delta\;{\overset{\_}{W}}_{PMSG}} = {\frac{1}{n}{\sum\limits_{i = 1}^{n}{\Delta{W_{PMSG}(i)}}}}$ with the discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) and its average value, obtaining the energy negative gradient by the following formula: $\eta = {- {\nabla\left\lbrack \frac{{{\Delta\;{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}}}{\left( {\frac{1}{n}{\sum\limits_{t = 1}^{n}\left( {{\Delta{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}} \right)^{2}}} \right)^{\frac{1}{2}}} \right\rbrack}}$ where ∇[⋅] is the symbol of gradient.
 4. The stability evaluation method of direct-drive wind turbine generator according to claim 3, assessing the system stability level according to η: when η is positive, the system is stable and the larger the value of η is, the higher the system stability level is; when η is zero, the system is critically stable; when η is negative, the system will go unstable, and as the absolute value of η increases, system stability level will drop.
 5. A stability evaluation system of direct-drive wind turbine generator, comprising: a Phasor Measurement Unit (PMU), a first processor, a second processor, a third processor and result output terminal, wherein the PMU is used to collect variations of voltage, current, active power and PLL (Phase Locked Loop) angle at a terminal of direct-drive wind turbine generator in one oscillation period $\frac{2\pi}{\omega};$ the first processor is used to calculate terminal energy of direct-drive wind turbine generator with a terminal energy model of direct-drive wind turbine generator based on the measured variations, wherein the terminal energy model of direct-drive wind turbine generator is: ${\Delta{W_{PMSG}(t)}} = {{\int_{t}^{t + \frac{2\pi}{\omega}}{\Delta P_{e}d\Delta\theta_{pll}}} + {\int_{t}^{t + \frac{2\pi}{\omega}}\left( {{\Delta i_{d}d\Delta u_{q}} - {\Delta i_{q}d\Delta u_{d}}} \right)}}$ where ΔP_(e) is the variation of direct-drive wind turbine generator output active power compared with a steady-state value, Δθ_(pll) is the variation of PLL angle, Δi_(d), Δi_(q), Δu_(q) and Δu_(d) are respectively the variation of d-axis and q-axis components of current and voltage at the terminal of direct-drive wind turbine generator, ω is the oscillation frequency of dq axis; the second processor is used to assess stability level of direct-drive wind turbine generator according to a value of energy negative gradient which is obtained according to discretized and average values of the terminal energy of direct-drive wind turbine generator, wherein the larger the value of the energy negative gradient is, the higher the system stability level is; the third processor is used to obtain an adjustment law of critical parameters affecting the stability of the system based on the energy negative gradient with stability operation points and critical parameters PLL and transmission lines; and the result output terminal is used to output the stability level of the system and the adjustment law of the critical parameters affecting system stability; wherein the third processor is operable for step S1: applying the PLL dynamic function to the energy negative gradient expression, achieving the energy negative gradient with stability operation points and PLL critical parameters: $\eta_{2} = {{- \left( {{I_{{sd}\; 0}\omega\frac{\cos\;\varphi_{c}}{2}} + {I_{{sq}\; 0}\omega\frac{\sin\;\varphi_{c}}{2}}} \right)}{\quad{{{\quad\left\lbrack {{{A_{-}(t)}\Omega_{u -}\alpha_{u}e^{\alpha_{u}t}} + {\frac{{dA}_{-}(t)}{dt}\Omega_{u -}e^{\alpha_{u}t}}} \right\rbrack }_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}} - {U_{{sd}\; 0}{\quad{\omega{\frac{\cos\;\left( {\varphi_{c} - \varphi_{RL}} \right)}{2}\left\lbrack {{{A_{-}(t)}\alpha_{i}\Omega_{i -}\ e^{\alpha_{i}t}} + {\frac{{dA}_{-}(t)}{dt}\Omega_{i -}\ e^{\alpha_{i}t}}} \right\rbrack}}}_{t_{0}}^{t_{0} + \frac{2\;\pi}{\omega}}} + {\left( {{I_{{sd}\; 0}\omega\frac{\cos\;\varphi_{c}}{2}} + {I_{{sq}\; 0}\omega\frac{\sin\;\varphi_{c}}{2}}} \right)\left\lbrack {{{A_{+}(t)}\Omega_{u +}\alpha_{u}e^{\alpha_{u}t}} + {\frac{d{A_{+}(t)}}{dt}\Omega_{u +}e^{\alpha_{u}t}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}{{+ U_{{sd}\; 0}}\omega{\frac{\cos\left( {\varphi_{c} - {\varphi\phi}_{RL}} \right)}{2}\left\lbrack {{{A_{+}(t)}\alpha_{i}\Omega_{i +}e^{\alpha_{i}t}} + {\frac{d{A_{+}(t)}}{dt}\Omega_{i +}e^{\alpha_{i}t}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}{{- 2}\;\omega\; K_{pll}^{2}\Omega_{u -}\Omega_{u +}e^{2\alpha_{u}t}\alpha_{u}\sin\;{\varphi_{c}\left\lbrack {{I_{sd0}{\sin\left( {\varphi_{u +} + \varphi_{u -}} \right)}} - {I_{sq0}{\cos\left( {\varphi_{u +} + \varphi_{u -}} \right)}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}}}}$ Where A⁻, A₊ respectively represent the amplitude of PLL angle variation Δθ_(pll) in sub-synchronous frequency and super-synchronous frequency; I_(sd0), I_(sq0), U_(sd0), U_(sq0) are the original value of dq axes components of current and voltage of the terminal of direct-drive wind turbine generator; α_(u) and α_(i) are the attenuation coefficients of oscillation components of voltage and current; Ω_(u−), Ω_(u+), Ω_(i−) and Ω_(i+) are respectively the amplitude coefficients of sub/super synchronous components of voltage and current; φ_(u−) and φ_(u+) are the initial phase angles of sub-synchronous and super-synchronous voltage; φ_(i−) and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents; Ω_(u+)=Ω_(i+)X₊, Ω_(u−)=Ω_(i−)X⁻, φ_(u−)=φ_(i−)+φ_(RL) and φ_(u+)=φ_(i+)+φ_(RL), where φ_(RL) is the phase angle of transmission lines; k_(pθ) and k_(iθ) are the proportion gain and integration gain of PLL; ${\varphi_{c} = {\arctan\frac{\omega k_{p\theta}}{k_{i\;\theta}}}},{{K_{pll} = {\frac{{k_{p\theta}s} + k_{i\;\theta}}{s^{2} - {U_{sq0}\left( {{k_{p\theta}s} + k_{i\;\theta}} \right)}}}_{s - {j\;\omega}}};}$ applying the dynamic function of system transmission lines to the energy negative gradient expression, obtaining the energy negative gradient of critical parameters of transmission lines: ${{\eta_{1} = {{- 2}\Omega_{i +}\Omega_{i -}e^{2\alpha_{i}t}\alpha_{i}X_{L}\omega\;{\sin\left( {\varphi_{i +} + \varphi_{i -}} \right)}}}}_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}} + {\quad{2\left\lbrack {{\Omega_{i -}^{2}e^{2\alpha_{i}t}\alpha_{i}\omega\; R} - {\Omega_{i +}\Omega_{i -}\ e^{2\alpha_{i}t}\alpha_{i}R\;\omega\;{\cos\left( {\varphi_{i +} + \varphi_{i -}} \right)}}} \right\rbrack}}_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}$ where Ω_(i−) and Ω_(i+) are the amplitude coefficients of current sub-synchronous component and super-synchronous component respectively; α_(i) is the attenuation coefficients of oscillation components of current; R and X_(L) are line resistance and the equivalent reactance of the transmission lines; φ¹⁻ and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents; and step S2: based on the constructed energy negative gradient with the critical parameters of stability operation points of direct-drive wind turbine generator, PLL and the transmission lines, obtaining influence pattern of these parameters on the energy negative gradient, wherein the adjustment law comprises: the stability level of system is improved by reducing I_(sd0) or increasing I_(sq0); the more stable the system is the larger R is and the smaller X_(L) is; the system stability is enhanced by reducing PLL parameters K_(p) and K_(i).
 6. The stability evaluation system of direct-drive wind turbine generator according to claim 5, the PMU equipped at the terminal of direct-drive wind turbine generator.
 7. The stability evaluation system of direct-drive wind turbine generator according to claim 5, where the second processor is used to assess the stability level of the system according to the value of energy negative gradient, comprising: S1: discretizing the terminal energy: ${\Delta{W_{PMSG}(i)}} = {\sum\limits_{t = t_{i}}^{t_{i} + \frac{2\pi}{\omega}}{\Delta{W_{PMSG}(t)}}}$ the value of terminal energy of direct-drive wind turbine generator at the time, t₁, . . . t_(n), i.e. ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) can be calculated; n is the number of data; S2: with discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)), calculating the energy negative gradient η, which comprises the following steps in detail: calculating the discretized average terminal energy of direct-drive wind turbine generator by the following formula: ${\Delta\;{\overset{\_}{W}}_{PMSG}} = {\frac{1}{n}{\sum\limits_{i = 1}^{n}{\Delta{W_{PMSG}(i)}}}}$ with the discretized terminal energy ΔW_(PMSG)(t₁), . . . , ΔW_(PMSG)(t_(n)) and its average value; obtaining the energy negative gradient can be based on the following formula; $\eta = {- {\nabla\left\lbrack \frac{{{\Delta\;{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}}}{\left( {\frac{1}{n}{\sum\limits_{t = 1}^{n}\left( {{\Delta{W_{PMSG}(t)}} - {\Delta\;{\overset{\_}{W}}_{PMSG}}} \right)^{2}}} \right)^{\frac{1}{2}}} \right\rbrack}}$ where ∇[⋅] is the symbol of gradient; and S3: assessing the stability level of system according to η: when η is positive, the system is stable and the larger the absolute value of η is, the higher the system stability level is; when η is zero, the system is critically stable; when η is negative, the system will go unstable, and as the absolute value of η increases, system stability level will drop.
 8. A non-transitory machine-readable storage medium comprising instructions that when executed cause a processor of a computing device to: measuring, by a Phasor Measurement Unit (PMU), variations of voltage, current, active power and PLL (Phase Locked Loop) angle at a terminal of direct-drive wind turbine generator in one oscillation period $\frac{2\pi}{\omega};$ based on the measured variations, calculating terminal energy of direct-drive wind turbine generator with a terminal energy model of direct-drive wind turbine generator, wherein the terminal energy model of direct-drive wind turbine generator is: ${\Delta{W_{PMSG}(t)}} = {{\int_{t}^{t + \frac{2\pi}{\omega}}{\Delta P_{e}d\Delta\theta_{pll}}} + {\int_{t}^{t + \frac{2\pi}{\omega}}\left( {{\Delta i_{d}d\Delta u_{q}} - {\Delta i_{q}d\Delta u_{d}}} \right)}}$ where ΔP_(e) is the variation of direct-drive wind turbine generator output active power compared with the steady-state value, Δθ_(pll) is the variation of PLL angle, Δi_(d), Δi_(q), Δu_(q) and Δu_(d) are respectively the variation of d-axis and q-axis components of current and voltage at the terminal of direct-drive wind turbine generator, ω is the oscillation frequency of dq axis; obtaining energy negative gradient of direct-drive wind turbine generator according to discretized and average values of the terminal energy of direct-drive wind turbine generator, and assessing system stability level according to a value of the energy negative gradient, wherein the larger the value of the energy negative gradient is, the higher the system stability level is; and obtaining an adjustment law of critical parameters affecting the stability of the system based on the energy negative gradient with stability operation points and critical parameters PLL and transmission lines, which comprises the following steps S1-S2: step S1: applying the PLL dynamic function to the energy negative gradient expression, achieving the energy negative gradient with stability operation points and PLL critical parameters: $\left. {\eta_{2} = {{- \left( {{I_{sd0}\omega\frac{\cos\varphi_{c}}{2}} + {I_{sq0}\omega}} \right)}\frac{\sin\varphi_{c}}{2}}} \right)$ ${\left\lbrack {{{A\_}(t)\Omega_{u -}\alpha_{u}e^{\alpha_{u}t}} + {\frac{{dA\_}(t)}{dt}\Omega_{u -}e^{\alpha_{u}t}}} \right\rbrack ❘_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}} -$ ${\left. {U_{sd0}\omega} \right){\frac{\cos\left( {\varphi_{c} - \varphi_{RL}} \right)}{2}\left\lbrack {{{A\_}(t)\alpha_{i}\Omega_{i -}e^{\alpha_{i}t}} + {\frac{{dA\_}(t)}{dt}\Omega_{i -}e^{\alpha_{i}t}}} \right\rbrack}❘_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}} +$ $\left. {{\left( {I_{sd0}\omega} \right)\frac{\cos\varphi_{c}}{2}} + {I_{sq0}\omega\frac{\sin\varphi_{c}}{2}}} \right)\left\lbrack {{{A_{+}(t)}\Omega_{u +}\alpha_{u}e^{\alpha_{u}t}} +} \right.$ $\left. {\frac{{dA}_{+}(t)}{dt}\Omega_{u +}e^{\alpha_{u}t}} \right\rbrack|_{t_{0}}^{{t_{0} + \frac{2\pi}{\omega} + {U_{sd0}\omega{\frac{\cos({\varphi_{c} - \varphi_{RL}})}{2}\lbrack{{{A_{+}(t)}\alpha_{i}\Omega_{i +}e^{\alpha_{i}t}} + {\frac{{dA}_{+}(t)}{dt}\Omega_{i +}e^{\alpha_{i}t}}}\rbrack}}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}} -$ ${2\omega K_{pll}^{2}\Omega_{u -}\Omega_{u +}e^{2\alpha_{u}t}\alpha_{u}\sin{\varphi_{c}\left\lbrack {{I_{sd0}{\sin\left( {\varphi_{u +} + \varphi_{u}} \right)}} - {I_{sq0}{\cos\left( {\varphi_{u +} + \varphi_{u -}} \right)}}} \right\rbrack}}|_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}$ Where A⁻, A₊ respectively represent the amplitude of PLL angle variation Δθ_(pll) in sub-synchronous frequency and super-synchronous frequency; I_(sd0), I_(sq0), U_(sd0), U_(sq0) are the original value of dq axes components of current and voltage of the terminal of direct-drive wind turbine generator; α_(u) and α_(i) are the attenuation coefficients of oscillation components of voltage and current; Ω_(u−), Ω_(u+), Ω_(i−) and Ω_(i+) are respectively the amplitude coefficients of sub/super synchronous components of voltage and current; φ_(u−) and φ_(u+) are the initial phase angles of sub-synchronous and super-synchronous voltage; φ_(i−) and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents; Ω_(u+)=Ω_(i+)X₊, Ω_(u−)=Ω_(i−)X⁻, φ_(u−)=φ_(i−)+φ_(RL) and φ_(u+)=φ_(i+)+φ_(RL), where φ_(RL) is the phase angle of transmission lines; k_(pθ) and k_(iθ) are the proportion gain and integration gain of PLL; ${\varphi_{c} = {\arctan\frac{\omega k_{p\theta}}{k_{i\theta}}}},{{K_{pll} = {❘\frac{{k_{p\theta}s} + k_{i\theta}}{s^{2} - {U_{sq0}\left( {{k_{p\theta}s} + k_{i\theta}} \right)}}❘}_{s = {j\omega}}};}$ applying the dynamic function of system transmission lines to the energy negative gradient expression, obtaining the energy negative gradient of critical parameters of transmission lines: $\eta_{1} = {{{- 2}\Omega_{i +}\Omega_{i -}e^{{2\alpha},t}\alpha_{i}X_{L}{{\omega sin}\left( {\varphi_{i +} + \varphi_{i -}} \right)}❘_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}} + {{2\left\lbrack {{\Omega_{i -}^{2}e^{{2\alpha},t}\alpha_{i}\omega R} - {\Omega_{i +}\Omega_{i -}\ e^{{2\alpha},t}\alpha_{i}R{{\omega cos}\left( {\varphi_{i +} + \varphi_{i -}} \right)}}} \right\rbrack}❘_{t_{0}}^{t_{0} + \frac{2\pi}{\omega}}}}$ where Ω_(i−) and Ω_(i+) are the amplitude coefficients of current sub-synchronous component and super-synchronous component respectively; α_(i) is the attenuation coefficients of oscillation components of current; R and X_(L) are line resistance and the equivalent reactance of the transmission lines; φ¹⁻ and φ_(i+) are the initial phase angles of sub-synchronous and super-synchronous currents; and step S2: based on the constructed energy negative gradient with the critical parameters of stability operation points of direct-drive wind turbine generator, PLL and the transmission lines, obtaining influence pattern of these parameters on the energy negative gradient; wherein the adjustment law comprises: the stability level of system is improved by reducing I_(sd0) or increasing I_(sq0); the more stable the system is the larger R is and the smaller X_(L) is; the system stability is enhanced by reducing PLL parameters K_(p) and K_(i). 